Pairing States of Superfluid ^He in Uniaxially Anisotropic Aerogel 
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Stable pairing states of superfluid "^He in aerogel are examined in the case with a global uniaxial 
anisotropy which may be created by applying an uniaxial stress to the aerogel. Due to such a 
global anisotropy, the stability region of an ABM pairing state becomes wider. In an uniaxially 
stretched aerogel, the pure polar pairing state with a horizontal line node is predicted to occur, as 
a 3d superfluid phase, over a measurable width just below the superfluid transition at Tc{P). A 
possible relevance of the present results to the case with no global anisotropy is also discussed. 

PACS numbers; 



In superconductivity and superfluidity, an intrinsic 
anisotropy has a profound effect on the resulting pair- 
ing state. For instance, one main origin of the d^-2_y-2 
pairing state in high Tc cuprates j3| is the four-fold 
symmetry of the Fermi surface resulting from the crys- 
talline anisotropy. Due to such an intrinsic and global 
anisotropy, a specific pairing symmetry with the highest 
temperature of Cooper instability is realized in the case 
of anisotropic superconductivity with no degeneracy be- 
tween different pairing states. In contrast, the bulk liquid 
■^He has an isotropic Fermi surface, and hence, there is a 
degeneracy in the transition point between different pair- 
ing states at least when the fluctuation is neglected [3;I3- 
Recently, possible pairing states of superfluid '^He in glob- 
ally isotropic aerogels were examined by comparing the 
free energy in the Ginzburg-LandaufGL) region between 
different states, and it was found Q that, after averag- 
ing over the quenched disorder brought by the aerogel 
structure, the pairing states to be realized are unaffected 
by the locally anisotropic scattering events due to the 
aerogel. However, situation may change if, as in the su- 
erconducting case mentioned above, a global anisotropy 
1^ is introduced externally in the aerogel structure. 

Here, we report on our results of theoretical phase 
diagrams of superfluid '^He in an aerogel with global 
anisotropy. Throughout this paper, we assume that such 
an anisotropy may be induced in quasiparticle scattering 
events by applying an uniaxial stress 0, Q to aerogels. 
Based on a conventional model Q of effects of aerogel on 
^He, the anisotropy can be incorporated in a momentum 
dependence of the random- averaged quasiparticle scat- 
tering amplitude. When the aerogel is uniaxially com- 
pressed, the 2d-like pairing state is favored, and the re- 
gion of the ABM pairing state becomes wider. In the 
stretched Id-like pairing state is favored at least 

near Tc, and actually the pure polar pairing state should 
be realized accompanied by a 2nd order transition to a 
deformed ABM state at a lower temperature. This is a 
rare situation in which a new pairing state is expected to 
occur as a 3d superfluid phase of '^He. This research was 
preliminarily reported elsewhere . 

Just as in Ref.Q, we start from the BCS hamiltonian 



with an impurity potential term 



rfV^Vi(i-)M(r)Va(r), 



(1) 



where u{v) denotes an impurity potential for quasiparti- 
cles brought by aerogel structures. As argued in Ref.j^j, 
the scattering amplitude |MkP in an aerogel, where Uk 
is the Fourier transform of u(r), should have a highly 
anisotropic and complicated momentum dependence, re- 
flecting local anisotropy in aerogels. However, as far as 
the aerogel is globally isotropic, the local anisot ropy is 
not reflected in averaged quantities such as |ukP, where 
the overbar denotes the average over the impurity con- 
figuration due to the aerogel. On the other hand, in a 
globally anisotropic aerogel, |ukP remains anisotropic. 
We invoke the model [g 



|ukP-A(l + <5„(k.z)^), 



(2) 



where 5u is a small parameter measuring the global 
anisotropy, k — k/Zcp, and z denotes the direction of 
a uniaxial deformation. Although the factor A is a func- 
tion of k^, such an isotropic k dependence induces no 
difference between various pairing states, and hereafter, 
A will be treated as a constant factor. Then, within the 
Born approximation, the quasiparticle Green's function 



Ge(p) = (ie-ep-SF 



satisfies 



^p-p' I 



le - 



Cp' ^p' 



(3) 



where e is the fermionic Matsubara frequency, and 
Jp, — J d^p' /{2ttY'. Below, we will focus on the situa- 
tions satisfying 2t:tT > 1, where r^^ = 27riV(0)(|ukp)k 
is the relaxation rate. Taking Ep to be purely imaginary 
and performing the p'-integral close to the Fermi surface, 
we obtain Sp = —ir/p sgne, where 



2r 



(l + 5„(p.z)2). 



(4) 



and Su ~ du{l + (5„/3)^^. The neglect of the real part, 
Re(I]p), of Sp is safely valid as far as 2tttT 3> 1: 
The p-dependence in Re(I]p) may be absorbed into an 



2 



anisotropy of density of states (DOS), while such a cor- 
rection to DOS is scaled by E^^ . Hence, the neglected 
correction to Re(I]) is of the order du/{EpT), which is 
smaller than the magnitude of particle-hole asymmetry 
Tc/E-p. In this manner, focusing on the imaginary part 
of E is justified. Then, Gslp) becomes 



Ge(p) = (i(|e| + ?7p)sgne - ^p)" 



(5) 



Note that, when Su > (< 0), the mean free path 
of normal quasiparticles running along the z-direction is 
shorter (longer). Thus, the case with i5„ > (< 0) cor- 
responds to the uniaxially compressed (stretched) case. 

Besides the self energy part, the vertex part in the 
particle-particle channel is also affected by the impurity 
scattering. If neglecting spatial variations of the pair-field 
A^^i, the bare vertex pi is replaced by Ti{e, p), where 

r.(e,p) = p.-f / r,(£,p')|G,(p')l'KwF- (6) 

Pl + {p- 



The solution of eq.® takes the form ri{e,p 
z)TgZi, where 



r. = -1 



1 - 2 
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k>l 



1 



l\2e 



(7) 



The above expressions will be used to derive a GL hamil- 
tonian per volume /iql in the anisotropic case. Up to 
0((5„), its quadratic term is expressed by 
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(8) 



where ipiz) is the di-gamma function, and v4^,i = 
|A(T)|a^,i with a^j,^i{a^j,^i)* — 1. The anisotropic term 
due to the vertex correction contributes to eq.© with 
the same sign as that due to Sp. 

The quartic term /iq^ is also derived in a similar man- 
ner and, up to 0((5„), takes the form 



' (4) 

'gl 



Each of the coefRcients /3i is the sum of a weak coupling 
contribution (3^*^ and a strong coupling one 5(3i. Regard- 
ing 5 Pi, their expressions with 5u — Q derived in Ref.0 
will be used hereafter. This approximation should not 
affect calculation results except for extremely large |(5„| 
values. The coefficients fif^ and /3|^^ are given by 
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up to 0(5„), where /3o(T) = 7C(3)A^(0)/(2407r2T2). 

As seen in eq.®, the inclusion of a global anisotropy 
induces a splitting of the Cooper instability point be- 
tween different pairing states. Since ip'-^^y) > (j/ > 0), 
an uniaxial compression with positive Su makes the insta- 
bility point of 2d-like pairing states with vanishing a^^z 
higher, implying that such a state must be realized just 
below Tc{P). Situation is similar to ^He thin films, and 
thus, this 2d-like state should be the ABM state. In the 
same manner, in the uniaxially stretched case with nega- 
tive Su, the Id-like polar pairing state with a^_i — d^Si^z 
tends to occur just below Tc. However, it is unclear at 
this stage whether or not the state just below may 
be a mixture of the ABM and polar pairing states Q so 
that the pure polar symmetry obtained at Tc crosses over 
upon cooling to the ABM one with no transition. 

It should be noted that, in the disordered case, ad- 
ditional terms are induced in the GL hamiltonian by 
the impurity scattering and its local anisotropy. When 
the global anisotropy is absent, a combination of one of 
such terms and the gradient term, which was not repre- 
sented in eq.(8), leads to destruction of superfluid long- 
ranged order (LRO) in the ABM state l|, T^. How- 
ever, the contribution to the free energy of the additional 
term is well described simply by incorporating a disorder- 
^i/jinduced shift of into the mean field condensation en- 
ergy i?c 0- In determining phase diagrams below, we 
have followed this finding j^] and, for brevity, have used 
the disorder- induced Tc-shift with 5„ = 0. The latter 
procedure does not affect our quantitative results unless 
all transitions between different pairing symmetries occur 
(9)in the close vicinity of Tc- 
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FIG. 1: Obtained P-T phase diagrams in an uniaxially- 
compressed case (a) with 5u = +0.04 and a stretched case 
(b) with (5„ = —0.07. For both figures, we have used 
(27rr)~^ = 0.137 (mK). The solid curves are transition lines 
in aerogel, while the thin dotted curves are those of bulk liq- 
uid. The ABM state in (a) is a genuine superfluid with LRO, 
while in (b) it is a superfluid glass 



We have numerically examined transitions between dif- 
ferent pairing states through by taking account of Tc- 
shifts of different origins mentioned above. Since fully 
including 18 real components of the pair field a^_i is 
cumbersome, nonvanishing four components, Rea^^^ and 
Imcj^y, were kept in calculations so that the familiar 
ABM. BW, planar, and polar states are taken into ac- 
count. For instance, in the uniaxially stretched case, 
the route through which only Rea^.z and Imoz^j, remain 
nonzero is found to be the most favorable upon cool- 
ing from the Cooper instability of the polar state (see 
Fig. 2). Typical examples of the resulting phase diagrams 
for 5„ > (compressed case) and i5„ < (stretched case) 
are given in Fig.l (a) and 1 (b), respectively, where the 
GL hamiltonian valid up to 0(^2) was used. In both of 
stretched and compressed cases, there is no polycritical 
pressure (PCP), i.e., there is a nonvanishing P-range of 
the ABM state even in low T limit. For the used r-value, 
there is no mean field stability region in P < 30 (bar) of 
the ABM state in the isotropic ((5„ = 0) case. It implies 
that the main origin inducing the ABM state in the fig- 
ures is the Tc-shift and not the strong coupling elfect. 
Further, we have verified that the planar pairing state 
cannot overcome the ABM one in free energy at any P 
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FIG. 2: Variation of a^,i at P = 22 (bar) around the 2nd 
order transition at Tpa between the polar and the deformed 
ABM states with |Reaz,z| 7^ |Imaz,j;|. The same material 
parameters were used as in Fig.l (b). A 1st order transition 
to the deformed BW state occurs at Tab- 



and T and is not realized as a pairing state in equilibrium. 

In the ABM state in the compressed case, the direction 
of gap nodes is pinned on average along z just like in thin 
films where the ABM or planar state has a wider stabil- 
ity region jj]. An origin of the remarkably wide ABM 
region in Fig.l (a) can be attributed to the similarity to 
the thin film case. Further, due to this pinning effect of 
1-vector, this ABM state has a true superfluid LRO in 
contrast to the quasi LRO in the isotropic case On 
the other hand, in the stretched case (or equivalently, the 
case compressed in the cylindrically symmetric manner) , 
the direction of gap nodes in the deformed ABM state 
(see Fig. 2) is spontaneously chosen within the x-y plane. 
Situation is similar to the bulk ^He in a uniform magnetic 
field which also favors the ABM state. Thus, the ABM 
stability region becomes wider even in the stretched case, 
although this state is a superfluid glass with no genuine 
superfiuid LRO Further, in both cases, the BW state 
with no gap nodes is deformed by the anisotropy (i.e, 

ax,x = ay,y 0.z,z) [13 • 

An intriguing result in the uniaxially-stretched case 
is the appearance of the pure polar pairing state, with 
a horizontal line of gap nodes in x-y plane, just below 
rc(P)-line. It appears that the temperature width over 
which the polar state is stable will be observable exper- 
imentally. Needless to say, this temperature width is, as 
well as that of the ABM state, extended with increasing 
\5u\- This polar state is not a mixture with other pairing 
states and, as in Fig. 2, shows a 2nd order transition to a 
deformed ABM state with point nodes in the x-y plane 
upon cooling. A Tc-shift dependent on the pairing states 
resulting from eq.(8) is essential to obtaining the pure 
polar state. Since the direction P of gap-maximum is 
pinned by z, this polar state has a true superfluid LRO. 

As a measure useful in detecting the polar state with 
the order parameter — Ad/^,, Pi, a pulsed NMR fre- 
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FIG. 3: Initial configuration to be taken in pulsed NMR ex- 
periments in the polar state. The angles 9 and /3 are the 
tipping angle in H - P plane and cos~^(P • H/ii"), respec- 
tively. 

quency shift in the polar state wiU be considered. Here, 
P points on averages to the stretched direction z. Fluc- 
tuations of P are assumed to be negligibly small. As 
usual, effects of the dipole energy on the spin dynamics 
can be examined in terms of Leggett's equations as 
far as the initial configuration is in equilibrium. In the 
pure polar state, the so-called dipole torque R^i 3 is 
given by 

RD = -^(d-P)(dxP), (11) 

5 

where qd is a dipole energy strength and hence, P _L 
d and H _L d are to be satisfied in equilibrium. Then, 
a frequency shift Acj of the free-induction signal due to 
R^) occuring by tipping the magnetization by an angle 9 
is given by 




(l + cose*)], (12) 



which depends on the angle (3 spanning P and H. 



Finally, we note that the present result may also be rel- 
evant to liquid ^He in globally isotropic aerogels ||3| if the 
correlation length of the local anisotropy 3 is much 
longer than the superfluid coherence length ^o- Since, as 
mentioned earlier, the free energy of each pairing state is 
roughly determined by the condensation energy Ec even 
at length scales of the order of ^q, the results induced by 
the global anisotropy, such as the wider ABM region 
and an occurrence of the polar pairing state near Tc, may 
be valid in the globally isotropic aerogel under the condi- 
tion 3> £,Q which may be satisfied at higher pressures. 
Then, a polar glass phase with spatially random P over 
long distances might be realized in the equal-spin pairing 
region near Tc. In contrast, at lower pressures with longer 
^o-values, the anisotropy-induced Tc-shift will be negligi- 
ble, and the approach ^3;] modelling the local anisotropy 
as a random field becomes appropriate. Then, a well- 
defined PCP is expected at a nonzero temperature. 

In conclusion, by introducing a global anisotropy in 
aerogels, the region of A-like phase should be extended if 
this phase has the ABM pairing state. A recent measure- 
ment has shown an extention of the A-like phase region, 
even on warming, due to an uniaxial compression 13]. In 
aerogels deformed via a uniaxial stretch or a cylindrically 
symmetric compression, an appearance of the pure polar 
pairing state near Tc and a wider region of ABM super- 
fluid glass Q are expected. We hope a measurement for 
the stretched case to be performed. 
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